Cosmography in /(T)-gravity 
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Being based on the only assumption that the universe is homogenous and isotropic on large 
scales, cosmography is an ideal tool to investigate the cosmic expansion history in a almost model- 
independent way. Fitting the data on the luminosity distance and Baryon Acoustic Oscillations 
allows to determine the confidence ranges for the cosmographic parameters hence giving some quan- 
titative constraints that a whatever theory has to fulfill. As an application, we consider here the 
case of teleparallel gravity (TEGR) also referred to as f(T) gravity. To this end, we first work 
out analytical expressions to express the present day values of f(T) derivatives as a function of 
the cosmographic parameters which hold under quite general and physically motivated conditions. 
We then use the constraints coming from cosmography to find out the confidence ranges for f(T) 
derivatives up to the fifth order and show how these can be used to check the viability of given 
TEGR models without the need to explicitly solve the second order dynamic equations. 

PACS numbers: 04.50.Kd, 98.80.-k 



I. INTRODUCTION 

Various cosmological observations, including the Type 
la Supernova (Sne la) the cosmic microwave back- 
ground radiation Q and the large scale structure [!, H| , 
ct al., have revealed that the universe is undergoing an 
accelerating expansion and it entered this accelerating 
phase only in the near past. This unexpected observed 
phenomenon poses one of the most puzzling problems in 
cosmology today. Usually, it is assumed that there exists 
dark energy (DE) in our universe, as an exotic energy 
component with negative pressure which dominates the 
universe filled with cold dark matter (CDM) and drives 
the universe to an accelerating expansion at recent times. 

The simplest and most appealing candidate for DE is 
the vacuum energy (cosmological constant, A) with a 
constant equation of state (EoS) parameter, — 1. This 
model is in general agreement with the current astro- 
nomical observations, but has difficulties to reconcile the 
small observational value of DE density to that coming 
from quantum field theories. This is called the cosmo- 
logical constant problem Q . Recently it was shown that 
ACDM model may also suffer from an age problem Q. 
It is thus natural to pursue alternative possibilities to 
explain the current acceleration of the universe. Observ- 
ing the small deviations for the EoS parameter from —1 
requires a description of the DE that allows the EoS to 
evolve across the phantom divide line —1 possibly mul- 



tiple times. The present data seem to slightly favor an 
evolving DE with the EoS parameter crossing —1 from 
above to below in the near past 0. One may take the 
observed accelerating expansion as a signal of the break- 
down of our understanding to the laws of gravitation and, 
thus, a modification of the gravity theory is needed. 
Over the past decade numerous DE models have been 
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proposed, such as quintessence 
[To| . tachyon [ll|, quintom 0, 

lyg in gas [r| and the generalized Chaplygin gas (GCG) 
14 1, the holographic DE [l5j . the new agegraphic DE 
16| , the Ricci DE 17[ and so on. 

In addition, the extensions to gravity by making the 
action a function of the spacetime curvature scalar R, 
f(R) [341H, or other curvature invariants 22 1, by cou- 
pling the Ricci scalar to a scalar field [23j], by introduc- 
ing a vector field contribution [24| 1 or by using properties 
of gravity in higher dimensional spacetimes have widely 
been investigated (25|. Among f(R) models, there are 
models that arc verified by all the observational and the- 
oretical constraints and exhibit universe acceleration and 
phantom crossing [26j-j30j. 

In a different approach, avoiding the curvature defined 
via the Levi-Civita connection, one could explore an al- 
ternative way and use the Weitzenbock connection that 
has no curvature but instead torsion. This has the prop- 
erty that the torsion is formed completely from products 
of first derivatives of the tetrad, with no second deriva- 
tives appearing in the torsion tensor. This approach 
was considered originally by Einstein in 1928 [3l|,|32|, as 
" Tclcparallclism" . It is closely related to standard Gen- 
eral Relativity, differing only in " boundary terms" involv- 
ing total derivatives in the action. The theory is naturally 
formulated by gauging external (spacetime) translation 
and underline the Weitzenbock spacetime characterized 
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by the metricity condition and by the vanishing of the 
curvature tensor. Translations are closely related to the 
group of general coordinate transformations which un- 
derlies General Relativity. The theory possesses a num- 
ber of attractive features both from the geometrical and 
physical viewpoints [33|-[48j]. 

Some models based on modification of the teleparallel 
equivalent of General Relativity (TEGR) are presented 
as an alternative to inflationary models without an in- 



a program has been particularly suited for the study of 
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flaton [4C . 14111 and DE models for the acceleration of the 



42 1-[48| where dark torsion (DT) is responsible 



universe 

for the observed acceleration of the universe, and the field 
equations are always 2nd order equations. This property 
makes these theories simpler than the dynamical equa- 
tions resulting in /(i?) gravity among other advantages. 
It has been shown that in theories of generalised TEGR, 
whose Lagrangians are algebraic functions of the usual 
teleparallel Lagrangian, the action and the field equa- 
tions are not invariant under local Lorentz transforma- 
tions . The authors also argue that the usual telepar- 
allel Lagrangian equivalent to General Relativity, is just 
a special case. 

It is worth noticing that all the above models such as 
dark energy, /(i?)-gravity and /(T)-gravity have shown 
to be, in broad sense, in agreement with observational 
data. As a consequence, unless higher precision probes 
of the expansion rate and the growth of structure will 
be available, these rival approaches could not be discrim- 
inated. This degeneration about the theoretical back- 
ground suggests that a more conservative approach to the 
problem of cosmic acceleration, relying on as less model 
dependent quantities as possible, is welcome. 

A possible solution could be to come back to the cos- 
mography rather than finding out solutions of the Fried- 
mann equations and testing them. Being only related 
to the derivatives of the scale factor, the cosmographic 
parameters make it possible to fit the data on the dis- 
tance - redshift relation without any a priori assumption 
on the underlying cosmological model. In this case, the 
only assumption is that the metric is the Robertson - 
Walker one. Almost a century after Hubble discovery of 
the expansion of the universe, we could now extend cos- 
mography beyond the search for the value of the Hub- 
ble constant. The SNela Hubble diagram extends up to 
z = 1.7 thus invoking the need for, at least, a fifth order 
Taylor expansion of the scale factor in order to give a 
reliable approximation of the distance - redshift relation. 
As a consequence, it could be, in principle, possible to 
estimate up to five cosmographic parameters, although 
the still too small dataset available does not allow to get 
a precise and realistic determination of all of them. 

Once these quantities have been determined, one could 
use them to put constraints on the models. In a sense, 
we are reverting the usual approach consisting in deriv- 
ing the cosmographic parameters as a sort of byprod- 
uct of an assumed theory. Here, we follow the other 
way around expressing the model characterizing quanti- 
ties as a function of the cosmographic parameters. Such 



fourth order theories of gravity, i.e. /(i?)-gravity [5 
As it is well known, the mathematical difficulties enter- 
ing the solution of fourth order field equations make it 
quite problematic to find out analytical expressions for 
the scale factor and hence predict the values of the cos- 
mographic parameters. But, nobody has still studied this 
procedure in f(T) gravity. A key role in both f(R) and 
/(T) scenarios is played by the choice of the f(R) or the 
f(T) function. Under quite general hypotheses, we will 
derive useful relations among the cosmographic param- 
eters and the present day value of f( n '(T) = d n f/dT n , 
with n = 0,...,5, whatever the analytic form of f(T) 
is. These relations will then allow to constrain the f(T) 
derivatives, provided model independent constraints on 
the cosmographic parameters are available. 

The layout of the paper is as follows. In Section II 
we will review the f(T) gravity. In Section III we will 
introduce the basic notions of the cosmographic param- 
eters. Section IV contains the main result of the paper 
demonstrating how the f(T) derivatives can be related 
to the cosmographic parameters, while, in Section V, we 
use these relations and previous constraints on the cos- 
mographic parameters to derive model independent esti- 
mates of the present day values of the f(T) derivatives. 
As a further application, Section VI shows how these lat- 
ter constraints can be used to observationally validate a 
given class of TEGR models without the need to solve the 
field equations. We summarize and conclude in Section 
VII. 



II. f(T) GRAVITY 

Teleparallelism uses as dynamical object a vierbein 
field ei(x M ), i — 0,1, 2, 3, which is an orthonormal basis 
for the tangent space at each point x^ of the manifold: 
ei.ej = r)ij, where rjij — diag(\, —1, —1, —1). Each vector 



ei can be described by its components e^, (i = 0, 1,2,3 
in a coordinate basis; i.e. = ef9 M . Notice that latin 
indexes refer to the tangent space, while greek indexes la- 
bel coordinates on the manifold. The metric tensor is ob- 
tained from the dual vierbein as g^ v {x) = nije l fi (x)el(x). 
Differing from General Relativity, which uses the torsion- 
less Levi-Civita connection, Teleparallelism uses the cur- 
vaturclcss Weitzcnbock connection, whose non-null tor- 
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(2.1) 



This tensor encompasses all the information about the 
gravitational field. The TEGR Lagrangian is built with 
the torsion (|2.1[) and its dynamical equations for the vier- 
bein imply the Einstein equations for the metric. The 
teleparallel Lagrangian is 
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where 



where 



S f r = \(K»% + <^T% - S;T 8 %) (2.3) 
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and K^ u p is the contorsion tensor 



K" u p = --(T^ p -T^ p -T p ^), 



(2.4) 



which equals the difference between Weitzenbock and 
Levi-Civita connections. 

In this work the gravitational field will be driven by 
a Lagrangian density that is a function of T. Thus the 
action reads 



I = 



1 



16ttG 



d A xe 



f(T), 



(2.5) 



where e = dei(e^) = y/—g. The case f(T) = T cor- 
responds to TEGR. If matter couples to the metric in 
the standard form then the variation of the action with 
respect to the vierbein leads to the equations [4(| 

e-%(eSinf(T) - e^T% x S p ^f{T) 

+Srd p (T)f\T) + \e»f{T) = 47rGe/T/,(2.6) 

where a prime denotes differentiation with respect to T, 
S^ v = ei p S p fll/ and is the matter energy-momentum 
tensor. 

We will assume a flat homogeneous and isotropic FRW 
universe, so 



diag(l, a(t), a(t), a(t)), 



(2.7) 



where a(t) is the cosmological scale factor. By using 
(EH), 4221), 4231) and 4H1) we obtain 
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(2.8) 



where H = - is the Hubble parameter. The substitution 
of the vierbein (|2.7[) in (|2 . 6[) for i = = v yields 

\2H 2 f'{T) + f(T) = WirGp. (2.9) 

Besides, the equation i = 1 = v is 

48H 2 f'(T)H - f'{T)[\2H 2 + 4H] - f(T) = !6nGp. 

(2.10) 

In Eqs. (|2.9[) and (|2.10p . p and p are the dark matter 
energy density and pressure, respectively. It can be easily 
derived that they accomplish the conservation equation 

p + 3H{p + P ) =0. (2.11) 

We can rewrite Eqs. (|2.9[) and (|2.10[) as the usual form 



H 



9 8ttG , 



2H + 3H 2 



8ttG 



(p + Pt) 



(2.12) 



(2.13) 



Pt 



16ttG 



[2Tf'(T) — f(T) — T/2], 



(2.14) 



1 



Pt 



16ttG 



[2H(4Tf"(T) + 2/'(T) - 1)] - p T . (2.15) 



are the torsion contributions to the energy density and 
pressure. Then, by using Eqs. (|2.14j) and (|2.15l) . we can 
define the effective torsion equation of state as 

_PT 1 , AH{ATf'{T) + 2f'{T)-l) 
U)T -y T =- l+ 4Tf>(T)-2f(T)-T ' (2 - 16) 

This could be, in principle, related to the observed accel- 
eration of the universe. 



III. COSMOGRAPHIC PARAMETERS 

Standard candles (such as SNela and, to a limited ex- 
tent, gamma ray bursts) are ideal tools in modern cos- 
mology since they make it possible to reconstruct the 
Hubble diagram, i.e. the redshift - distance relation up 
to high redshift values. It is then customary to assume 
a parameterized model (such as the concordance ACDM 
one, or any other kind of dark energy scenario) and con- 
trasting it against the data to check its viability and con- 
straints its characterizing parameters. As it is clear, such 
an approach is model dependent so that some doubts al- 
ways remain on the validity of the constraints on derived 
quantities as the present day values of the deceleration 
parameter and the age of the universe. In order to over- 
come such a problem, one may resort to cosmography, i.e. 
expanding the scale factor in Taylor series with respect to 
the cosmic time. Such an expansions leads to a distance 
- redshift relation which only relies on the assumption 
of the Robertson -Walker metric thus being fully model 
independent since it does not depend on the particular 
form of the solution of cosmic equations. To this aim, it 
is convenient to introduce the following functions: 



H 

q 
j 

s 

I 



1 da 

a dt 

a dt 2 
a dt 3 
adt 4 

i.*±H-* 

a dt 5 



a H~ 



-4 



(3.1) 
(3.2) 
(3.3) 
(3.4) 
(3-5) 



which are usually referred to as the Hubble, decelera- 
tion, jerk, snap and lerk parameters, respectively. Their 
present day values (which we will denote with a subscript 
0) may be used to characterize the evolutionary status of 
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the universe. For instance, qo < denotes an acceler- 
ated expansion, while j$ allows to discriminate among 
different accelerating models. 

It is then a matter of algebra to demonstrate the fol- 
lowing useful relations: 

H = -H 2 (l + q), (3.6) 
H = H 3 (j + 3q + 2), (3.7) 
H = H 4 [s - 4j - 3q(q + 4) - 6], (3.8) 

H (iv) = H 5 [l-5s + 10{q + 2)j + 30(q + 2)q + 24], (3.9) 

where a dot denotes derivative with respect to the cosmic 
time t and = d 4 H/dt 4 . Eq_s.§^ - make it 

possible to relate the derivative of the Hubble parameter 
to the other cosmographic parameters. 

IV. /(T) DERIVATIVES VS COSMOGRAPHY 

Rather than choosing a parameterized expression for 
f(T) and then numerically solving modified Friedmann 
equations for given values of the boundary conditions, we 
try to relate the present day values of its derivatives to 
the cosmographic parameters (qo,jo,so,lo) so that con- 
straining them in a model independent way gives us a 



hint for what kind of f(T) model could be able to fit the 
observed Hubble diagram. 

As a preliminary step, it is worth considering Eq. (|2.8[1 . 
Differentiating with respect to t, we easily get the follow- 
ing relations: 

t = -12HH, (4.1) 

f = -12[H 2 + HH], (4.2) 

f = -12[3HH + HH], (4.3) 

T (iv) = ^i2[ZH 2 + AH'ii + HH {m) ] (4.4) 

The modified Friedmann Eqs. (|2.9[) and (|2.10|) can be 
rewritten as 

H 2 1 = Y^ f - ) [Tn m + f(T)] (4.5) 

and 

H = JJJJfjl™ m - 4Hff"(T)] (4.6) 

where dot denotes derivative with respect to the cosmic 
time t and fl m represents the dimcnsionlcss matter den- 
sity parameter. However, in order to enter other cosmo- 
graphic parameters we have to differentiate from Eq. (|4.6j) 
three more times. We thus get: 



and 



H = 4g "r (T) \HT - T(3H 2 + 2H)] - J^li^f + HT)f"(T) + HT 2 f"'(T)}, (4.7) 

H = ; r^?,,^ [T{9H 4 + 6H 2 H + AH 2 ) - HT(3H + 6H 2 ) + H(HT - 2HT)] 
4H z f'(T) 

- JjJTojT, [HHf'(T) + {2H 2 f + 3HHT + 4HHT + H 2 f)f"(T) + H 2 f 3 f (m) {T) 

+ Hf(4Hf + 3Hf)f'"(T)}, (4.8) 

H {m] = - \T(10HHH + 12H 3 H - 27 H 6 - 12H 2 H 2 - 8H 3 - 2H 2 H) + H 3 f 

AH 3 f'(T) L v ; 

+ H 2 f(9HH + 27 H 3 - 5H) - 3H 2 f(3H 2 + H) + 7HH 2 f] 

■ [(3HHH + H 2 H + HH 2 )f'(T) + H 2 f 2 (7HT + 6HT)f {lv) (T) 



H 2 f'{T) 1 

+ {AH 2 HT + 2H 3 f + 7H 2 H'f + 10HH 2 T + 1H 2 HT + 11HHHT + H 3 T (lv) )f" \T) 
+ H(lOH 2 f 2 + 7HHT 2 + 21HHTT + 3H 2 f 2 + AH 2 ff)f'"{T) 

+ H 3 f 4 f {v \T)], (4.9) 
with /M(T) = d i f(T)/dT i and f {v) (T) = d 5 f(T)/dT 5 . 



Since the last five equations have to hold along the full evolutionary history of the universe, they naively hold 
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also at the present day. As a consequence, we may eval- 
uate them in t = t thus easily obtaining: 



Hl = 



-1 



12/'(T ) 



[T n m0 + /(T )] 



Hn = 



4/'(t ; 



[T Q m0 - AH Q T f"(T )] 



(4.10) 



(4.11) 



and so on for the next three ones. 

So, we have five equations, i.e. (|4.10|) . (|4.11j) and (|4~7| - 
(|4.9[) evaluated at the present day. We call these, "final 
equations" which will turn out to be useful in the follow- 
ing. But, one further relation is needed in order to close 
the system and determine the six unknown quantities 
/(To), /'(T ), /"(T ), /'"(To), /(™)(T ) and /W(T ). 
This may be easily obtained by noticing that, inserting 
back the physical units, equation (|2.9|) reads: 



H 



8ttG 
6/'(T) 



/cm 



16ttG 7 ' 



(4.12) 



Let us now suppose that /(T) may be well approxi- 
mated by its fifth order Taylor expansion in T — To, i.e. 
we set: 

/(T) = /(To) + /'(To)(T-To) + i/"(To)(T-T ) 2 
+ lf"'(T )(T - T ) 3 + i/(»)(T )(T - To) 4 
+ ^/^(ToXr-To) 5 . (4.15) 

In such an approximation, it is f( n '(T) — d n f/dT n = 
for n > 6. 

Evaluating Eqs.(|4.ip - and ((2~51) at the present 

time and using Eqs. (|3.6[) - (|3.9|) . one gets: 



Ti 



-6H, 



o • 



To = 12fl£(l+<R>), 



(4.16) 



(4.17) 



which clearly shows that, in /(T) gravity, the Newtonian 
gravitational constant Gpj has to be replaced by an ef- 
fective (time dependent) coupling G e ff. However, the 
present day value of Newtonian gravitational constant 
has to be recovered, and then: 



G eff (z = 0)=G N ^ f'(T ) = l 



(4.13) 



which means the recovery of TEGR. In other words, the 
choice /(T) = T gives rise to T = —6H 2 (see Eq. 
and then Eq. (|4.12[) reduces to the standard 



H 



8nG 



(4.14) 



To = -12ftf[«fo(g + 5) + j + 3], (4.18) 
To = -12i? 5 [ So - jo(3g + 7) - 3<z (4<Zo + 9) - 12], (4.19) 



T M = -12H^[l - s (4?o + 9) + io(3jo + 44 9o + 48) 



3g (4gg + 39 9o + 56) 



(4.20) 



After inserting all of these into the "final equations" , we 
can solve them under the constraint (|4. 1 3[) with respect 
to the present day values of f(T) and its derivatives up 
to the fifth order. After some algebra, one ends up with 
the desired result: 



/(To) 
67T 2 



mO 



f'(T ) = 1, 



f"(T ) -3O m0 1 
(eiT 2 )- 1 4(1 + go) 2' 



(4.21) 
(4.22) 
(4.23) 



/'"(To) _ -3fi m o(3(7 2 + 6go+Jo + 2) , 3 



(6i? 2 )- 2 



(1 + QoY 



+ 4' 



(4.24) 



f {iv) (T ) 

(e^o 2 )" 3 



mO 



16(1 + go) 5 



[*o(l + 9o ) + io(6g§ + 17g + 3j + 5) 



3g (5gg + 20gg + 29g + 16) + 9] 
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(4.25) 



6 



(6ff 2 )- 4 



32(1 + g ) 

„2 



? [io(l + qaY + so(10<z5 + 43gg + 46g + 13) + 10jb*o(l + So) 

32) 



+ 5j 2 (6g^ + 22q + 3jo + 7) 
+ 3gb(35g§ + 210g£ + 518gg 



■io(45^ + 225gg + 412gg 



666gg 



448go + 150) + 60] + — 



219g 4 
105 

16 ' 



(4.26) 



r 



Egs. (|4.2ip - (|4.26p make it possible to estimate the 
present day values of f(T) and its first five derivatives 
as function of the Hubble constant H and the cosmo- 
graphic parameters (qo,jo, So> lo) provided a value for the 
matter density parameter f2 m o is given. 



mographic parameters for ACDM model as simplest case. 
This is a minimal approach but it is useful to probe the 
self-consistency of the model. The cosmographic param- 
eters for the ACDM model read 



A. The ACDM model 

In order to get a first hint on the possible values of 
f(T) and its derivatives we have to reproduce the cos- 



J 



s = -(^) 2 (i-o m0 



H 



2a 3 ' 



= (^) 3 (l-O m0 + ^f) 3 / 2 , 



H 



2a 3 +"m0+ 2a 3 2a 6 J' 



(4.27) 
(4.28) 
(4.29) 



ti 



fi 2 

1 s ro0 



osi m0 



350 2 

° Ji 'mO \ / 

2a 6 jl 



1 - a 



a 



(4.30) 



which, evaluated at the present time, give 

3 

<?0 = — 1 + TT^mO, 



./() 



so — 1 — 2^ m0 ' 



27 

Zo = 1 + 3Sl m o H — — 



?n0 ' 



(4.31) 
(4.32) 
(4.33) 

(4.34) 



Inserting the previous equations in Eqs. (|4.23|) - (|4.26l) . 
we obtain 

f"(To) = /'"(To) = fW(T ) = f^(To) = 0, (4.35) 

and in the absence of these terms f(T) reduces to f(T) ~ 
T — 2 A. This is consistent with what we expected for the 
ACDM model and can be assumed as a consistency check. 



V. OBSERVATIONAL CONSTRAINTS 



In order to constrain the model by cosmography, i.e. 
to estimate the function /(T) through its own value and 
that of its derivatives at the present time, we need to 
constrain obscrvationally the cosmographic parameters 
by using appropriate distance indicators. Moreover, we 
must take care that the expansion of the distance re- 
lated quantities in terms of (qo, jo, sq, h) closely follows 
the exact expressions over the range probed by the data 
used. Taking SNela and a fiducial ACDM model as a 
test case, one has to check that the approximated lumi- 
nosity distance 1 deviates from the ACDM one less than 
the measurement uncertainties up to z ~ 1.5 to avoid in- 



See 15011 for the analytical expression. 



7 






-0.4 -0.3 -0.2 -0.1 0.0 0.1 0.2 



-0.3 -0.2 -0.1 00 0.1 0.2 



-0.3 -0.2 -0.1 00 0.1 0.2 



FIG. 1: Isolikelihood (68, 95 and 99% CL) contours for the /< quantities. The fuzzyness is due to numerical artifacts. 



x xbf (x) x m ed 68% CL 95% CL 

h 0.718 0.706 0.706 (0.693, 0.719) (0.679, 0.731) 

go -0.64 -0.44 -0.43 (-0.60, -0.30) (-0.71, -0.26) 

jo 1.02 -0.04 -0.15 (-0.88, -0.90) (-1.07, 1.40) 

s -0.39 0.18 0.02 (-0.57, 1.07) (-1.04, 1.78) 
l 4.05 4.64 4.54 (2.99, 6.48) (1.78, 8.69) 



x xbf (x) x med 68% CL 



95% CL 



fo -1.742 -1.733 -1.733 (-1.743, -1.723) (-1.751, -1.712) 

/a -0.033 0.113 0.147 (0.007, 0.208) (-0.153, 0.226) 

f 3 -0.092 0.530 0.815 (0.172, 0.921) (-1.483, 1.033) 

U 0.294 -0.955 1.061 (0.193, 2.306) (-18.307, 3.603) 

,/s 8.690 -68.893 6.371 (2.956, 11.014) (-370.966, 31.004) 



TABLE I: Constraints on the cosmographic parameters. 
Columns are as follows : 1. parameter id; 2. best fit; 3., 
4. mean and median from the marginalized likelihood; 5., 6. 
68 and 95% confidence ranges. 



TABLE II: Constraints on the /j values from the Markov 
Chain for the cosmographic parameters. Columns order is 
the same as in Tabled] 



troducing any systematic bias. Since we are interested in 
constraining (<7o,io> s o ; 'o)i we Wlu expand the luminosity 
distance Dl up to the fifth order in z which indeed allows 
us to track the ACDM expression with an error less than 
1% over the full redshift range. We have checked that 
this is the case also for the angular diameter distance 
Da = Dl(z)/(1 + z) 2 and the Hubble parameter H(z) 
which, however, we expand only up to the fourth order 
to avoid introducing a further cosmographic parameter. 

In order to constrain the parameters (h,qo,jo,So,lo), 
Bouhmadi - Lopez et al. [ll| have used the Union2 SNela 
dataset |HH and the BAO data from the analysis of the 
SDSS seventh release (55| adding a prior on h from the re- 
cent determination of the Hubble constant by the SHOES 
team [56j]. We update here their analysis adding the 
measurement of H{z) obtained in 1 5 211 from the age of 
passively evolving galaxies and in [53[ from the radial 
BAO. Exploring the five dimensional parameter space 
with a Markov Chain Monte Carlo method, we obtained 
the constraints summarized in TableU in agreement with 
previous results in literature [57T4591 ] . Note that, because 



of the degeneracies among the five cosmographic param- 
eters, the best fit values can also be different from the 
median ones, which is indeed what happens here. This 
is, however, not a shortcoming of the fitting analysis, but 
a consequence of the Bayesian approach giving more im- 
portance to sampling the marginalized parameters distri- 
butions rather than to looking for the best fit accordance 
within a given model and the available dataset 2 . In par- 
ticular, here, the best fit values are quite close to those 
predicted for the ACDM model (for instance, jo = 1 for 
a A dominated universe), while the median ones allow 
for significant deviations (with the ACDM values being, 
however, within the 95% confidence ranges). 

In order to translate the constraints on the cosmo- 



2 Qualitatively, one can say that the best fit value of, e.g., go is l es s 
important than the median one since the best fit go is the correct 
one only if the other parameters also take their best fit values, 
while the median one is more reliable since it describes the full 
distribution whatever are the values of the other parameters. 
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graphic parameters on similar constraints on the present 
day values of f(T) and its derivatives, we should just use 
Eqs. (|4.21j) - (|4.26j) evaluating them along the final coad- 
ded and thinned chain of the cosmographic parameters 
and then looking at the corresponding histograms. To 
this end, however, we should set the value of fi m o which 
is not constrained by the fitting analysis described before. 
To overcome this difficulty, we rely on the WMAP7 deter- 
mination of the physical matter density uj rn = £l m oh 2 = 
0.1329 and, for each value of h along the chain, we fix 
fl m o = uj m /h 2 having neglected the error on to m since 
it is subdominant with respect to the one on h. Note 
that the adopted estimate of w m comes from the fit to 
the CMBR anisotropy spectrum and mainly depends on 
the early universe physics only. Since it is reasonable to 
expect that GR is recovered in this limit, we can safely 
assume the validity of this result whichever f(T) model 
is considered. Defining for shortness 



A. Dependence on the expansion order 

Although we have checked that our fifth order expan- 
sion closely matches the exact luminosity and angular di- 
ameter distances and the Hubble parameter within less 
than 1%, it is worth noting that a decent approximation 
is also obtained if we stop the expansion to the third or 
fourth order. Cutting the expansion to order three (four) 
means that we can only constrain cosmographic param- 
eters up to the jerk jo (the snap So) and hence work out 
confidence ranges for the f(T) derivatives up to the third 
(fourth) order. It is nevertheless worth exploring how the 
constraints depend on the order of the expansion. To this 
end, we fit the same dataset as above with both the third 
and fourth order expansion of the involved quantities and 
then use the corresponding Markov Chains to estimate 
confidence limits on (/o, f%, fs). From the third order fit, 
we get (median and 68 and 95% CL) : 



fn = f {n \T )/(6H 2 )-^ 



we finally get the constraints summarized in TablclPIl and 
shown in Fig. [1] where the degeneracy between some cou- 
ples of parameters is shown as an example. Note that, 
as best fit value, we mean the one obtained by fixing the 
cosmographic parameters to the best fit values. However, 
because of the degeneracies among (qo,ja, sq,Iq) and the 
nonlinear behavior of the relations with f n , it is possible 
that the best fit /„ are quite different from their median 
values which is indeed the case (in particular, for f§). 
Note also that the confidence ranges become larger as 
the order n of the derivative increases. This is indeed 
an expected result since the higher is n, the larger is 
the number of cosmographic parameters involved so that 
the weakness of the constraints on the higher order cos- 
mographic parameters and the degeneracies among them 
makes the constraints on f n weaker and weaker as n gets 
larger. From a different point of view, such a behavior 
simply reflects the naive expectation that one has to go 
to deeper redshifts to probe the exact functional shape 
of f(T) and hence put severe constraints on the value of 
its high order derivatives. As a further remark, we note 
that the constraints on (/3, / 4 , fa) are strongly asymmet- 
ric with a long tail extending towards negative values 
causing a large offset between the mean and the median. 
This is actually a consequence of the term (l+go) _Q > with 
a = (3,5,7) for (/3,/4,/s) respectively, which enters as 
a common factor in Eqs. (|4.24[) - (|4.26|) . As q comes close 
to -1, these term becomes increasingly large so thus mak- 
ing /„ explode. It is, however, wort noting that values 
of qo close to -1 are indeed quite unlikely (although still 
allowed by our fit to a limited dataset) so that only the 
95% confidence ranges are affected. 



f _ -I 741 +0.009 +0.017 
J0— 1- '^-o. 008 -0.016 



n nn^+ 054 +0.098 

U.UUO_ a0 gg _ .154 



t — n nQ7+ - 303 +0-475 



552 



while the fourth order fit gives : 



fa = -1.733+ 011 + - 021 



h = 0.043 



-0.009 -0.019 



+0.061 +0.113 
-0.076 -0.195 



t — fl 4QQ+0.266 +0.439 
I .13 - U.4dy_ 441 _!.53 6 

Comparing the value of fi for the different fits (including 
the fifth order one in Table HI)) allows us to draw some 
interesting lessons. First, although the median values 
are different, the confidence ranges are well overlapped 
thus indicating that the order of the expansion should 
not have any statistically meaningful impact on the con- 
straints. However further and accurate studies have to 
be performed in order to confirm this statement. On the 
other hand, increasing the order of the expansion shifts 
away from the ACDM one (i.e., fi = for i > 1). This 
is, actually, a subtle effect of the degeneracy among the 
cosmographic parameters. Indeed, increasing the order 
n of the expansion adds further parameters to the fit 
thus allowing for much more combinations of the cosmo- 
graphic parameters able to fit well the same data. As a 
consequence, the constraints on go w iH become weaker al- 
lowing for models with q$ closer to —1 and hence (/2, /s) 
values far away from the fiducial ACDM ones. We, nev- 
ertheless, recommend the user to refer to the results in 
TableHIl since the fifth order expansion provides a bet- 
ter approximation to the underlying expansion history 
so that the fit is less affected by any bias due to any 
error in the approximation. 
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B. Deviations from basic assumptions 



Imposing Eq. (|4.2ip and f'(T ) = 1 gives: 



The constraints discussed above have been obtained 
under two basic underlying assumptions. First, we have 
set /'(T ) = 1 in order to recover an effective gravita- 
tional constant which matches the Newton one today. 
Actually, although reasonable, there are no compelling 
reasons why the Newton constant which is measured in 
laboratory experiments is the same as the cosmological 
one. As such, it is worth wondering how our results 
would change should we allow for deviations from the 
Gyv = G cosmo assumption. 

On the other hand, we have used the WMAP7 con- 
straints on the physical matter density ujm to infer 
the present day matter density parameter and then use 
Eqs.gjn]) - to constrain (/ , fa, fa, k, fa) from the 

cosmographic parameters. Some recent works [60l ] have, 
however, investigate the evolution of perturbations in 
/(T) theories finding out remarkable differences with re- 
spect to the standard GR. As a consequence, one can 
not exclude the possibility to recover a correct growth 
of structure even if /(T)) does not reduce to GR in the 
early universe. Should this be the case, the use of the 
WMAP7 ujm value is incorrect. 

Taking care of these possible effects is actually quite 
easy Indeed, some algebra shows that Eqs. (14.211) - 14.261 
can all be recast as : 



2 - n m0 - [1 + (H m0 - 2)6} In Tp 
l + (<y-l)lnT 



(f!mO — 1)^0 ' 

1 + In To 



(6.2) 



(6.3) 



so that we can express fa for i = (2, 3, 4, 5) as function of <5 
only. We then proceed as follows. For each fa value of the 
sample obtained above from the cosmographic parame- 
ters analysis, we solve fa (5) = fa- Since this equation has 
two roots, we store them and then compute (fa, f 4, fa) 
for both values thus obtaining an histogram for the model 
prediction of these quantities. The median and 68% and 
95% confidence ranges read : 



f _ n 9Qfi+ - 272 +0-599 
J3 - -U-^9O_ . lls _ .149 



f n SQ1 +0-330 +0.424 

J4 — U.OJ1_ 799 _i.797 

f — _o c;«a+3.143 +7.176 
^ J5 — 0.dOO_ 1 274 _ x 633 

choosing the lowest 5 solution and 



, Vn g0, JO Jo, So) , 

(1 + q ) a " 

with V n (qo, jo, lo, s o) a polynomial function of its argu- 
ments, a„ = (0, 1, 3, 5, 7) for n = (0, 2, 3, 4, 5), n n a con- 
stant depending on n and we have set /'(To) = 1 + e. 
Using this simple formula allows to immediately scales 
our constraints to different values of the ujm and e pro- 
vided one has a theoretical or observational estimate of 
these quantities. 



VI. COSMOGRAPHY VS f(T) MODELS 

Up to now, we have never assumed any functional 
shape for f(T) so that the constraints in Table I indeed 
holds for the full class of TEGR theories provided one can 
approximate f(T) by its fifth order Taylor series over the 
redshift range probed by the data. Such a result can also 
be read in a different way. Given a f(T) model, its char- 
acterizing parameters must be chosen in such a way that 
the constraints in Table I are satisfied. This considera- 
tion offers an interesting route to check the viability of a 
given f(T) model without the need of explicitly solving 
the field equations and fitting the data. 

As an example, let us assume the following model [6l| : 



f — a 77Q+0.193 +0.415 
J3 — O.I m_ oms _n.131 



-3.120 



-0.018 
-0.024 



-0.032 
-0.050 



-31.033 



+0.371 +0.525 
-0.811 -1.810 



for the larger solution. Since the 95% CL in TablellTl 
are quite large because of the impact of qo, wc will use 
only the 68% confidence ranges which we compare the 
above constraints to. For the lower 5 solutions, both fa 
and fa are smaller than the 68% CL from cosmographic 
parameters, while the range for fa has a marginal overlap. 
On the other hand, choosing the largest 8 solution leads 
to [fz, fa, f$) values that fully disagree with the model 
independent constraints. We therefore argue that the 
mode l (16. ip is disfavored by the observational data. 
In [61( , another model was also proposed : 



f(T) =aT + l3T n 



(6.4) 



where, imposing as before the constraints on /(To) and 
/'(T ), one easily gets: 



(2 - n m0 )n - 1 
n-1 



(6.5) 



f(T) = aT + PT s \nT 



(6.1) 



(flmO — 1)T 

n-1 



l-n 



(6.6) 
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We then solve fa(n) = fa and estimate the theoretically 
expected values for the other derivatives obtaining : 



n oct;+0.272 +0.599 
" u - zod -0.116 -0.148 



t — n Q/1 1 +0.331 +0.424 

f _ _o oi 7+3.156 +7.215 
U5- °- OJ -' -1.282 -1.641 



These values are still in disagreement with the constraints 
in Tablelrfl hence making us argue against this model 
too. Actually, some caution is needed in this case. If 
we set a = 1 and \n\ small enough, Eq. (|6.4j) predicts 
an expansion rate which can be made arbitrarily close 
to the ACDM one. Indeed, if we use the best fit value 
of the cosmographic parameters, we find n = —0.011 
and quite small values for {fajf^fa) as expected for A 
term. Actually, the disagreement with the constraints in 
TablclITl may be due to a failure of one of the underly- 
ing assumptions in the derivation of Eqs. (|4.21l) - (|4.26|) . 
Indeed, these relations have been obtained by Taylor ex- 
panding f(T) to the fifth order thus implicitly assuming 
that the higher order terms are subdominant. Depending 
on the value of n, however, this assumption can fail for 
the model (|6.4I) so that the constraints on /„ should not 
be considered reliable. 



VII. CONCLUSIONS 

Cosmography offers a valid tool to investigate cosmic 
expansion in a model independent way. The constraints 
on the cosmographic parameters (go, jo, sq, Iq) obtained 
by fitting to SNela Hubble diagram and BAO data are 
fully general relying on the only assumption that the 
universe is homogenous and isotropic on large scales. 
As such, any given cosmological model should predict 
(<7o> jo, sq, Iq) values which are in agreement with these 
constraints. Such a premise makes it clear why studying 
the cosmography of a given theory can offer a valuable 
help to check its viability as an explanation of the ob- 
served cosmic speed up. 

Motivated by these considerations, we have discussed 
the cosmography of TEGR theories obtaining the expres- 
sion of f(T) and its derivatives as a function of the matter 
density parameter f2 m o and the cosmographic parame- 
ters (h,q ,j ,s Jo)- It is worth stressing that the re- 
lations thus found hold for all TEGR models provided 
they can be well approximated by their fifth order Tay- 
lor expansion, at least over the redshift range probed by 
the data used to constrain the cosmography. A key role 
has been played by the assumption /'(To) = 1, mean- 
ing that the effective gravitational constant equals the 
Newton one at redshift z = 0. Although this is a quite 
reasonable assumption, it nevertheless rely on the under- 
lying identification of the cosmological G cosmo with the 
local one Gn- Should /'(To) 7^ 1, one could rc-derive 



our results, but the price to pay is to lower the order of 
the expansion of one degree. Alternatively, one can leave 
fa = /'(T ) as a free parameter and check which is the 
impact on (fa, / 3 , fa, fa). Not surprisingly, if fa ~ 1 + e 
with |e| << 1, one could still use our relations to con- 
strain (fa, fa, fi, fa) making a systematic error which is 
by far smaller than the statistical uncertainties unless 
unreasonable large values |e| (> 0.1) are adopted. 

The above relations allow to transform the con- 
straints on (n m Q,h) and the cosmographic parameters 
(go, jo, «o, '0) into similar ones for the (fa, / 3 , fa, fa) quan- 
tities. Coming out from a model independent approach 
as cosmography, these constraints have to be fulfilled by 
any TEGR model. As such, we can investigate a priori 
(i.e., without solving the field equation) of a given f(T) 
theory by simply comparing the theoretically predicted 
(/2, / 3 , fa, fa) with the observed ones. As an application, 
we have considered here two particular classes showing 
that, although they can in principle give rise to an ac- 
celerated expansion, they are both unable to predict the 
observationally motivated {fa, fa, fa, fa) values so that 
they can be rejected. Expanding on this idea, one could 
also reverse the approach and build up a class of theories 
that fits the above constraints from the beginning and 
then investigate which is the expansion history at higher 
z and the growth of structure. 

However, we have to stress again that the method pro- 
posed here has some shortcomings: i) the truncation of 
the Taylor expansion at some predefined order can be 
problematic due to the dropped higher-order terms. Such 
terms could reveal important and so an arbitrary trun- 
cation could be dangerous; ii) one would face difficulty 
when trying to retain as high an order as possible, be- 
cause the high-order terms are increasingly weakly con- 
strained by the limited data and could have large errors 
(for example, the constraints given in Table II and those 
in Sec. VA could not give /(T)-curves which agree very 
well with each other) . In this case, the second shortcom- 
ing is made worse by the first one. 

To be more specific, the new constraints on the 
(ho,qo,jo,so,lo) parameters make sense, as well as their 
translation to the fa (i = 0, 2, • • • ,5) parameters, consid- 
ered the fact that the higher i is, the bigger and more 
uncertain is fa. This is a weakness of the method used 
here, namely the higher-order terms in the Taylor expan- 
sion, which could be important in the overall behavior of 
f(T), are more difficult to predict accurately. It is worth 
noticing that smaller deviations of (ho, 90, jo, s o, '0) from 
the corresponding ACDM values could cause very big de- 
viations of fas from their corresponding ACDM values 
which are identically zero. 

Another comment is in order for the results in Sec. VA. 
The fits have been performed with different orders of the 
Taylor expansion. We can see that the fa parameter 
perfectly agrees with the value in Table II, but not fa or 
fa. The reason is that we have used a different number of 
parameters to fit the same curve. The assumption is that 
the order of expansion has no statistically meaningful 
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impact on the constraints, but such a statement should 
be confirmed by further studies using more complete data 
sets. On the other hand, it is clear that the increasing 
order of expansion shifts away from the ACDM fiducial 
values. 

Although these drawbacks, the approach is interesting 
and might be made more accurate as soon as more data, 
especially those coming from highcr-redshift surveys, will 
be taken into account. 

As a final comment, it is worth noticing how the re- 
newed interest in old dated cosmography has now opened 
the way to an alternative and yet powerful method to in- 
vestigate, on the same ground, both dark energy models 



and modified gravity theories, such as f(R) and TEGR 
models. After so many years, however, we are no more 
interested in finding only two numbers, namely (h,qo), 
but rather we now need a fifth order expansion, hence 
five quantities, to constrain not only the evolution of the 
universe, but also the underlying theory of gravity. 
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